Low energy graphene edge termination via small diameter nanotube formation 
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We demonstrate that free graphene sheet edges can curl back on themselves, reconstructing as 
nanotubes. This results in lower formation energies than any other non-functionalised edge struc- 
ture reported to date in the literature. We determine the critical tube size and formation barrier 
and compare with density functional simulations of other edge terminations including a new recon- 
structed Klein edge. Simulated high resolution electron microscopy images show why such rolled 
edges may be difficult to detect. Rolled zigzag edges serve as metallic conduction channels, separated 
from the neighbouring bulk graphene by a chain of insulating sp^-carbon atoms, and introduce Van 
Hove singularities into the graphene density of states. 



The atomic structure of graphene and graphene edges, 
is a subject of great interest, in particular stimulated 
by new aberration corrected atomic resolution electron 
microscopy studies |l|, |2|. Edge structure can define 
the chemical and electronic properties of graphene rib- 
bons [3], yet there is no consensus about the most stable 
free edge structure. Unterminated edges, consisting in 
a line of atoms with dangling bonds, are inherently un- 
stable and subject to chemical functionalisation in ambi- 
ent. However under vacuum conditions, such as in elec- 
tron microscope columns, unterminated edges can be ob- 
served. Simply cutting through the graphene lattice re- 
sults in the most studied edge structures: armchair or 
zigzag, or a combination of the two. The first way that 
free edges can stabilise themselves is via rehybridisation 
of the carbon atoms. This occurs spontaneously for the 
armchair edge giving a sequence of double bonds along 
the edge. The zigzag edge has been shown to be also 
met ast able and undergo a 5-7 reconstruction y, which 
has been recently experimentally confirmed [5] . High res- 
olution transmission electron microscopy (HRTEM) sug- 
gests that alternative edge structures may be also com- 
mon, including the theoretically less stable zigzag and 
Klein edges [2, 6-11]. Other reconstructed edges, loosely 
based on the Haeckelite structures [12| have also been 
proposed [1|. Besides reconstruction, HRTEM [2|, |6|- 
lll Il3l Il4| has shown that free edges can fold back on 
themselves resulting in no longer a graphene monolayer 
(so-called grafold [15]). The energetic cost of bending the 
layer is partially compensated by Van der Waals interac- 
tions in the stacked region. 

In this Letter we extend the range of edge types, in- 
troducing a third type of stabilised edge whereby edges 
are rolled back on themselves and rebonded into the 
graphene sheet. This results in a graphene monolayer 
with nanotube at the edge site, and eliminates all dan- 



gling bonds by sp^-like rehybridisation of the carbon 
atoms along the rebonding line. 

We perform spin polarized density functional calcula- 
tions under the local density approximation using the 
AIMPRO code [16]. Graphene edges are modeled using 
ribbons with width ~50A in orthogonal supercells large 
enough to avoid interaction between neighbouring cells. 
Tubular rolled edges were created symmetrically on both 
of the two ends of the graphene ribbon. Edge formation 
energies per unit length (eV/A) are determined via 
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where E^ot is the total internal energy of a system with n 
carbon atoms, jic is the energy of a carbon atom in a per- 
fect graphene sheet and L is the ribbon edge length. Re- 
action barriers were calculated using the climbing nudged 
elastic band technique [17], with all atoms allowed to 
move. More details are given in the supplementary ma- 
terials. 

We begin by considering formation energies for zigzag, 
armchair, 5-7 reconstructed zigzag, Klein [18] and recon- 
structed Klein edges (Fig. [T] and supplementary mate- 
rials). In all cases no out of plane distortion (e.g. edge 
rippling) was found after relaxation. In agreement with 
results of Ref. [j|, we obtain the armchair and 5-7 recon- 
structed zigzag edges to be the most energetically stable. 
The least stable is the Klein edge with unsaturated car- 
bon atoms, however if symmetry is allowed to break, it 
spontaneously reconstructs by rebonding in pairs, gain- 
ing ~0.7 eV/A (see supplementary materials). Segments 
of unreconstructed Klein edge were recently observed by 
annular dark field microscopy p, ll9| but the observed 
Klein edges might have residual hydrogen not detectable 
by annular dark field microscopy. 

A new type of termination leading to the stabiliza- 
tion of the edge can be achieved by taking an unrecon- 
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Figure 1. (Color online) Formation energy (eV/A) for dif- 
ferent graphene edges (dashed lines indicated on the figure: 
Klein, reconstructed Klein, zigzag, armchair and 5-7 recon- 
structed zigzag), free standing zigzag and armchair tubes, 
armchair-, armchair- like-, and zigzag-tube terminated edges. 




Figure 2. (Color online) Upper panel. Structure of (4,4) arm- 
chair (a), (4,4) armchair- like (b) and (8,0) zigzag (c) nan- 
otube terminated graphene sheet, sp^-like coordinated carbon 
atoms are marked in yellow. Bottom panel. Corresponding 
simulated HRTEM images. 



structed free edge, folding it back on itself and bond- 
ing the edge dangling bonds to a line of basal graphene 
atoms. A rolled zigzag edge can bond into the graphene 
plane in two configurations: either above what would 
be zigzag edge atoms or Klein edge atoms. Rolled arm- 
chair edges can only bond to equivalent armchair edge 
atoms. Drawing on nanotube nomenclature we refer to 
these new edges as armchair-nanotube terminated (Fig. 
[2^), armchair-like nanotube terminated (Fig. Wp)^ and 
zigzag-nanotube terminated (Fig. [2t) respectively. In 
all cases the line of carbon atoms bridging the tube and 
graphene layer adopts an sp^-like hybridization with av- 
erage bond lengths (~ 1.50 A) and angles (~ 108°) close 
to those of diamond. Locally, the structure is similar 
to the core of the zigzag prismatic dislocation in AA 
graphite [20]. The sp^-like bonding allows the tube to 
localise strain, resulting in a droplet-shaped cross-section 
(Fig. 0. 

Fig. [1] presents formation energies for different types 
of free and tube terminated edges as a function of the 



Hamada index of the tube. For comparison we plot the 
values for free-standing armchair and zigzag tubes. Con- 
sistent with an earlier study [21|, free standing small 
tubes with diameters below 4 A (i.e. (3,3) and (5,0) and 
below) are unstable compared to a flat graphene sheet. 
This means that it is thermodynamically preferable to 
split these tubes open, even with unfunctionalised edges. 
Indeed experimentally these small radii nanotubes have 
not been observed on their own but they can exist as 
inner tubes in large multiwalled nanotubes [22|. We see 
here however that such small nanotubes are more stable 
when formed on graphene ribbon edges, through locali- 
sation of the nanotube strain along the sp^-coordinated 
tetrahedral bonding line, e.g. a (3,3) armchair tube and 
all zigzag tubes up to (7,0). Comparing to free zigzag, 
edges become more stable when rolled in tubes above 
(8,8). Rolled armchair edges have a lower energy than 
free edges when forming nanotubes larger than (14,0). 
For the largest presented tube terminated edges, forma- 
tion energies are lower than any of the previously pro- 
posed free edge configurations. 

If we extrapolate to larger diameter tubes we might ex- 
pect that they collapse due to Van der Waals interactions 
between walls to a dog-bone cross-section [23|. However 
the droplet cross-section induced by the line of sp^-like 
carbon atoms naturally induces a "local collapse", and 
this pinched region extends further as the tube diameter 
increases. Thus for large diameters, rolled edges converge 
to a classical folded edge which then terminates some dis- 
tance from the actual edge via a line of sp^-like bonds. 
We note that in no case is the combination of a free tube 
and graphene edg e more stable than the nanotube ter- 
minated edge [2^ i.e. there will be no thermodynamic 
driving force for nanotube production from rolled edges. 

We next examined the barrier to roll up the edges, de- 
termining the barrier between the free zigzag edge and an 
armchair-nanotube terminated edge as a function of the 
final tube size. For a (4,4) nanotube terminated edge, we 
find an NEB barrier for ribbon rolling of about 2 eV/A, 
with a corresponding unrolling barrier of 0.9 eV/A. For 
a (8,8) nanotube terminated edge barriers change to 1.3 
eV/A to roll versus 1.6 eV/A for derolling. Increas- 
ing further the tube diameter gives tube reconstructed 
edges more stable than free edges and this induces a 
barrier height inversion. In particular for small tubes 
these barriers are quite high, due to the high curvature 
that has to be induced. However these represent max- 
imum barriers and experimental barriers will likely be 
much smaller, since our calculations assume concerted 
simultaneous bonding/debonding along the entire edge 
length. Similar to dislocation motion, which proceeds via 
the propagation of kinks along the dislocation line rather 
than a concerted single-step motion, rolled edges will pre- 
sumably roll/unroll initially at a single point which will 
then propagate along the nanotube length. 

Given their low formation energies compared to free 
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Figure 3. Calculated density of states for (a) (8,8) armchair- 
and (b) (8,0) zigzag-nanotube terminated graphene edge. The 
Van Hove singularities are marked by arrows. 



zigzag edges, we can ask the question why such tube ter- 
minated edges have not yet been reported in the litera- 
ture. One reason is that these edges require long-range 
order, whereas other edge structures can vary over the or- 
der of single unit cells. Thus we may expect such edges 
to be more common in well defined periodic ribbon edges 
such as after splitting of large mult i- walled nanotubes. 
We note also that free standing graphene can be obtained 
through wet etching after epitaxial growth on metal sub- 
strate ^]. We expect that during the etching process 
rolled edges might appear as in analogous synthesis mech- 
anisms used for the production of inorganic nanotubular 
materials [26]. Furthermore, rolled edges could be diffi- 
cult to discriminate by transmission electron microscopy. 
Simulated HRTEM images (Fig. [2]) are very similar to 
those of free standing edges, the primary difference be- 
ing minor variations in the image contrast. Additionally 
these edges may display characteristic nanotube modes 
detectable using spatially resolved resonant Raman. 

We next examine the effect of rolled edges on the elec- 
tronic properties of the graphene. We present in Fig. 
[3] the electronic density of states for the (8,8) and (8,0) 
tube terminated edges. These show an interesting combi- 
nation of the graphene and nanotube behaviour. In spite 
of different electronic character of free standing (8,8) and 
(8,0) tubes, i.e. metallic and semiconducting, in both 
cases the composite system has a non-zero density of 
states at the Fermi level. For a zigzag tube termination 
the background density of states around the Fermi level 
rises smoothly, reflecting the graphene density of states, 
overlaid on which there is a series of Van Hove singulari- 
ties characteristic of a nanotube. For the armchair- tube 
terminated edge, i.e. a rolled zigzag graphene edge, there 
is a sharp peak at the Fermi level similar to that seen for 
flat unterminated zigzag edges [sl, Ij . Such Fermi level 
peaks can lead to magnetic instability, and indeed our 
spin unrestricted calculations shows a ferromagnetic con- 



figuration to be the lowest energy state (slightly lowering 
of the edge energy of the system by only 0.02 eV/A), al- 
though this must be treated with caution given the use 
of the local density approximation. 
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Figure 4. (Color online) (a) Band structure for (8,8) armchair- 
nanotube terminated edges, and (b) distribution of highest 
occupied state at the F-point (up) and X-point (down). 

The peak consists of two degenerate states which can 
be seen in the associated band structure (Fig. |4^) as in- 
volving a mixing of several bands. Plotting these states 
at the F and X point (Fig. IIJd) shows that the Fermi level 
spike is localised mainly on the row of graphene atoms 
next to the sp^-carbon atoms, mirroring the edge state 
seen in fiat zigzag terminated graphene [3|, whereas at 
the X-point the highest occupied state comes from dis- 
persive edge states located in the nanotube segment along 
the junction. This convergence of three zigzag edges at 
the line of sp^-bonded carbon atoms suggests possible 
interesting magnetic behaviour under an applied field. 

We note that the zigzag metallic edge state is preserved 
in this configuration, and unlike the unterminated sim- 
ple zigzag edge it will also be partially protected from 
environmental attack, since all neighbouring atoms are 
fully coordinated. Notably we might expect it to be sta- 
ble in air. These states also suggest intriguing transport 
behaviour, with possible conduction channels both along 
the edge of the graphene and in the edge states in the 
nanotube segment. 

In summary we have shown that by rolling an unter- 
minated graphene edge it is possible to create nanotube- 
terminated edges where the sheet edge rebonds back into 
the graphene plane. We determine the critical tube size 
and formation barriers and compare with density func- 
tional simulations of other edge terminations including a 
new reconstructed Klein edge. We find that the proposed 
tube terminated edges are more stable than any other 



non-functionalized edge structure, due to the replace- 
ment of dangling bonds with sp^-like hybridised carbon 
atoms. Rolled zigzag edges serve as metallic conduction 
channels, separated from the neighbouring bulk graphene 
by a chain of insulating sp^-carbon atoms, and intro- 
duce Van Hove singularities into the graphene density of 
states. They may provide a way to stabilise and protect 
from chemical attack the disperse Fermi level state seen 
along metallic zigzag edges. Similar edge rolling effects 
might also appear in other layered materials [27] such as 
boron nitride monolayers where orbital rehybridization 
can occur. 

This work has been carried out within the NANOSIM- 
GRAPHENE project nANR-09- NANO-016-01 funded 
by the French National Agency ( ANR) in the frame of its 
2009 programme in Nanosciences, Nanotechnologies and 
Nanosystems (P3N2009). 



* V . ivanovskaya@gmail . com] 
^ chris.ewels@cnrs-imn.fr 

[I] X. Jia, J. Campos-Delgado, M. Terrenes, V. Meunier, and 
M. S. Dresselhaus, Nanoscale 3, 86 (2011). 

[2] C. Girit, J. Meyer, R. Erni, M. Rossell, C. Kisielowski, 
L. Yang, C. Park, M. Crommie, M. Cohen, S. Louie, et aL, 
Science 323, 1705 (2009). 

[3] T. Enoki, Y. Kobayashi, and K. Fukui, Int. Rev. Phys. 
Chem. 26, 609 (2007); J. Kunstmann, C. Ozdogan, 

A. Quandt, and H. Fehske, Phys. Rev. B 83, 045414 
(2011); T. Wassmann, A. Seitsonen, A. Saitta, M. Lazzeri, 
and F. Mauri, Phys. Rev. Lett. 101, 96402 (2008). 

[4] P. Koskinen, S. Malola, and H. Hakkinen, Phys. Rev. 

Lett. 101, 115502 (2008). 
[5] P. Koskinen, S. Malola, and H. Hakkinen, Phys. Rev. B 

80, 73401 (2009); A. Chuvilin, J. Meyer, G. Algar a- Siller, 

and U. Kaiser, New J. Phys. 11, 083019 (2009). 
[6] Z. Liu, K. Suenaga, P. J. F. Harris, and S. lijima, Phys. 

Rev. Lett. 102, 015501 (2009). 
[7] J. H. Warner, F. Schaffel, M. H. Riimmeli, and 

B. Biichner, Chem. Mat. 21, 2418 (2009). 

[8] J. Y. Huang, F. Ding, B. L Yakobson, P. Lu, L. Qi, and 

J. Li, Proc. Natl. Acad. Sci. USA 106, 10103 (2009). 
[9] M. Cass, U. Bangert, A. Bleloch, P. Wang, R. Nair, and 

A. Geim, Nat. Nanotechnol. 3, 676 (2008). 

[10] J. Meyer, A. Geim, M. Katsnelson, K. Novoselov, 
T. Booth, and S. Roth, Nature 446, 60 (2007). 

[II] S. Rotkin and Y. Gogotsi, Mat. Res. Innov. 5, 191 (2002). 
[12] H. Terrones, M. Terrones, E. Hernandez, N. Grobert, 

J. Charlier, and P. Ajayan, Phys. Rev. Lett. 84, 1716 
(2000). 
[13] J. Warner, M. H. Riimmeli, A. Bachmatiuk, and 

B. Biichner, Nanotechnology 21, 325702 (2010). 

[14] J. C. Meyer, A. K. Geim, M. L Katsnelson, K. S. 

Novoselov, D. Obergfeh, S. Roth, C. Girit, and A. Zettl, 

Solid State Commun. 143, 101 (2007). 
[15] K. Kim, Z. Lee, B. Malone, K. T. Chan, B. Aleman, 

W. Regan, W. Gannett, M. F. Crommie, M. L. Cohen, 

and A. Zettl, Phys. Rev. B (2011), in Press. 
[16] R. Jones and P. Briddon, Semicond. Semimetals 51A, 



287 (1998); M. Rayson and P. Briddon, Comput. Phys. 
Commun. 178, 1 28 (2008); M. J. Rayson and P. R. Brid- 
don, Phy s. Rev. B 80, 205104 (2009) , 

[17] G. Mills and H. J0nnson, Phys. Rev. Lett. 72, 1124 
(1994). 

[18] D. Klein, Chem. Phys. Lett. 217, 261 (1994). 

[19] K. Suenaga and M. Koshino, Nature 468, 1088 (2010). 

[20] L Suarez-Martinez, G. Savini, G. Haffenden, J. Campan- 
era, and M. Heggie, Phys. Status Solidi C 4, 2958 (2007). 

[21] S. Sawada and N. Hamada, Solid State Commun. 83, 917 
(1992). 

[22] L. Qin, X. Zhao, K. Hirahara, Y. Miyamoto, Y. Ando, 
and S. lijima, Nature 408, 50 (2000); L. Guan, K. Sue- 
naga, and S. lijima, Nano Lett. 8, 459 (2008). 

[23] N. Chopra, L. Benedict, V. Crespi, M. Cohen, S. Louie, 
and A. Zettl, Nature 377, 135 (1995). 

[24] L. Chernozatonskii, E. Sheka, and A. Artyukh, JETP 
Lett. , 352 (2009). 

[25] S. Shivaraman, R. Barton, X. Yu, J. Alden, L. Her- 
man, M. Chandrashekhar, J. Park, P. McEuen, J. Parpia, 
H. Craighead, et a/., Nano Lett. 9, 3100 (2009). 

[26] O. Schmidt and K. Eberl, Nature 410, 168 (2001). 

[27] J. Coleman, M. Lotya, A. G'Neih, S. Bergin, P. King, 
U. Khan, K. Young, A. Gaucher, S. De, R. Smith, et a/., 
Science 331, 568 (2011). 



Supplementary materials. 

Computational Details Structures were geometrically 
optimized using a (1 x 1 x 12) shifted k-point grid gener- 
ated using the Monkhorst-Pack formalism. Hartwigsen, 
Goedecker and Hutter relativistic pseudopotentials were 
usedpy]. Atom-centered Gaussian basis functions are used 
to construct the many-electron wave function with angu- 
lar momenta up to 1=2, i.e. 22 independent functions 
per C atom. We confirmed that all structures were en- 
ergetically converged with respect to ribbon width by 
extending some of the calculations by another hexagon 
row, finding formation energy changes corresponding to 
the addition of bulk graphene atoms. Electronic level oc- 
cupation was obtained using a Fermi occupation function 
with kT = 0.04 eV. Density of states calculations used a 
denser (1 x 1 x 100) k-point grid, with a 0.03 eV Gaussian 



broadening applied. 

HRTEM image simulations High resolution transmis- 
sion electron microscopy (HRTEM) images were sim- 
ulated using a full dynamical calculation by multislice 
method as implemented in the SimulaTEM program [2|. 
Formation energies and HREM image simulation for dif- 
ferent edge types are presented in Fig. [5l A focal series 
for a (8,8) armchair rolled tube is presented in Fig. [6l 
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Figure 5. Calculated formation energies for different edge structures and corresponding simulated HRTEM images. 
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Figure 6. HRTEM simulated focal series images for (8,8) arm- 
chair tube terminated edge. 



